Primordial beryllium as a big bang calorimeter 
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Many models of new physics including variants of supersymmetry predict metastable long-lived 
particles that can decay during or after primordial nucleosynthesis, releasing significant amounts of 
non-thermal energy. The hadronic energy injection in these decays leads to the formation of 9 Be via 
the chain of non-equilibrium transformations: Energy h — > T, 3 He — > 6 He, 6 Li — > 9 Be. We calculate 
the efficiency of this transformation and show that if the injection happens at cosmic times of a few 
hours, the release of 0(10 MeV) per baryon can be sufficient for obtaining a sizable 9 Be abundance. 
The absence of a plateau-structure in the 9 Be/H abundance down to a 0(lO~ 14 ) level allows one to 
use beryllium as a robust constraint on new physics models with decaying or annihilating particles. 



The rare light elements lithium, beryllium, and boron 
play a pivotal role in observational cosmology. Their 
abundance pattern is key in our understanding of stellar 
structure and galactic evolution. It may also probe the 
very beginning of nuclear chemistry, known as big bang 
nucleosynthesis (BBN). Standard BBN (SBBN) theory 
has become what could be called a parameter-free theory. 
With the cosmic microwave background determination 
of the baryon asymmetry, r) b ~ 6.23 x lCT 10 1], SBBN 
solely relies on Standard Model physics embedded in a 
Friedman-Robertson- Walker Universe. An overall con- 
cordance between the SBBN predictions and the obser- 
vationally inferred primordial abundances of D and He 
has put BBN on firm footing and made its framework an 
invaluable toolbox for testing models of new physics (see, 
e.g. 0). Whereas Li, Be, and B are mainly produced in 
galactic cosmic rays [3,111]; the only stable isotope to reach 
out to a primordial fraction at the ppb level is 7 Li. The 
SBBN yields of Be and B are negligible [j| and, to date, 
the only known scenarios which may produce Be in inter- 
esting amounts are inhomogeneous nucleosynthesis, see 
e.g. [6j, and the catalysis of nuclear reactions by electro- 
magnetically or strongly interacting relic particles ■ 

The purpose of this paper is to show that non- 
equilibrium processes induced by the decay or annihi- 
lation of a relic particle species X after BBN can 
lead to the production of primordial 9 Be/H at the 
0(1O -14 ) level and above, which makes these scenarios 
testable via atmospheric measurements of 9 Be in ex- 
tremely metal-deficient stars (MDSs). We identify the 
following prospective path to 9 Be: 



Energy^ T, 3 He -> 6 He, 6 Li -)• 9 Be, 



(1) 



where Energy^ refers to the injection of energy via 
hadronic channels. The chain of transformations is ini- 
tiated by the spallation of 4 He due to "primary" nu- 
cleons and followed by repeated reactions on ambient 
a-particles. It is well known that such non-thermal 
cascades — truncated at A = 6 — can be an efficient source 
of "secondary" 6 Li 1(| as it bypasses the extremely 
inefficient SBBN channel D(cv,7) 6 Li. Here we show that 



the exoergic reaction 

6 He + 4 He -> 9 Be 



Q = 0.60 MeV 



(2) 



is particularly capable of going one step further and 
bridging the A = 8 divide. As we shall argue, 9 Be has the 
potential of becoming one of the "calorimeters" of choice 
in detecting non-thermal hadronic processes during this 
early epoch of our Universe. We also investigate a similar 
"tertiary" production of 10 B. 

Observations The elements Li, Be, and B are de- 
tected via atomic resonance lines in the photospheres 
of MDSs. The constancy of 7 Li/H at lowest metallici- 
ties [Fe/H] < —1.5 is known as the "Spite plateau" [l2j ]. 
and is interpreted as representing the primordial lithium 
component. Whereas numerous determinations exist — 
with 7 Li/H in the range (1 — 2.5) x 10~ 10 — the observa- 
tional status of the isotopic ratio 6 Li/ 7 Li remains doubt- 
ful. Though a plateau value of 6 Li/ 7 Li ~ 0.05 has been 
claimed from detections in a number of halo stars [Dij . 
this claim has also been challenged. Since the 6 Li absorp- 
tion line is not resolved spectroscopically from the one of 
7 Li, convective motions may well mimic its presence [T3 |. 
Indeed, a re-analysis [H| has diminished the number of 
such (3tr) detections to a mere number of two stars. 

In contrast to Li, Be and B show no indication of a 
plateau structure. Instead, they exhibit a scaling with 
metallicity, which points towards a pure cosmic ray ori- 
gin. The lowest observationally inferred B abundance 
then sets a limit on its primordial value, (B/H)p < 
5 x 10~ 12 16]. The determination of beryllium is espe- 
cially robust since 9 Be is the only stable isotope. FigureQ] 
shows a sample of 9 Be detections (taken from [TtJ) as a 
function of metallicity. The lowest point in black corre- 
sponds to the upper limit obtained recently in [3] and 
which we are going to use as an indicative upper limit on 
primordial 9 Bc: 



( 9 Be/H) p < 10 



-14 



(3) 



If instead only the detected values of 9 Be are used, a nom- 
inal 3cr upper limit of 4.7 x 10 -14 can be deduced. In any 
case, 9 Be is by far the rarest and the most constrained 
among the light elements. 
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FIG. 1. Beryllium detections (squares) and fit of their Be-Fe 
trend taken from [l7l |. The triangle shows the upper limit 
from 0]; A(X) = log 10 (X/H) + 12. 



In addition to being firmly detected, there lies an- 
other virtue of Be and B. 6 Li is fragile and may well 
have been destroyed in MDSs. Such a destruction is 
indeed the essence of the "astrophysical solution" to 
the 7 Li problem, which posits a uniform stellar deple- 
tion of 7 Li by a factor 3 to 5 from its SBBN prediction 
7 Li/H = (5.24t° ^)xl0^ 10 19]. Any mechanism achiev- 
ing this requires that the atmospheric material is trans- 
ported deep enough into the hotter stellar interiors where 
7 Li is destroyed. If all three elements are subjected to the 
proton-burning reactions at the same temperature, one 
obtains 



7 Li/ 7 Li p = ( 6 U/ 6 U p ) a = ( 9 Be/ 9 Be p 



(4) 



with the exponents given by a = (cr^Li+pv) / ( <76 Li+p w ) 
and b = (o~7 Li+p v) / {<79 Bc+p v) (which makes the esti- 
mate exponentially sensitive to temperature.) Evaluat- 
ing the cross sections for 6 Li(p, a) 3 He, 7 Li(p, a) 4 He, and 
9 Be(p, a) 6 Li as well as 9 Be(p, D) 8 Be at a fiducial temper- 
ature of T = 10 6 K one finds a ~ 0.01 and b ~ 3 x 10 5 . 
Irrespective of a rather crude nature of (U), this esti- 
mate illustrates one point: had 7 Li been depleted by 
an 0(l)-factor, 6 Li could have been suppressed substan- 
tially, while 9 Be is to remain unaffected [20j. Therefore, 
if stellar depletion is indeed a prevalent factor in the re- 
duction of 7 Li, the use of 9 Be over 6 Li for constraining 
non-standard BBN scenarios may be preferred. Table U 
summarizes the observational status of the rare light el- 
ements 6 Li, 9 Be, and B = 10 B + 11 B, where we make a 
generous allowance for the possible stellar depletion of 
primordial 6 Li. 

Conversion probability The key quantity that deter- 
mines the efficiency of (JlJ and ultimately the strength of 
9 Be constraints is Pbc, that we define as the probability 
for an energetic A = 3 nucleus injected into the BBN 
plasma to produce 9 Be. Clearly, the chain ([1]) can be 
broken up into its sub-processes, 



Pbc - (P3->6) E in * (Pq^e™ * { P 9^9} 



TABLE I. Rare light elements and (a) their SBBN abundance 
prediction (b) observationally inferred upper limit on the pri- 
mordial abundance and (c) the temperature below which their 
BBN destruction by protons effectively stops; B = 10 B + n B. 



element 


SBBN 


primordial limit 


T f /keV 


6 Li/H 


~ lO" 14 


< few x 10" 10 


8 


9 Be/H 


< 10~ 18 


< 10" 14 


10 


B/H 


< 10" 15 


< 5 x 10~ 12 


20 



where each individual probability Pa 1 ^a 2 is averaged 
over the distribution of At injection energies 25™ . It is 
easy to see that Pb b has to be a functional of the initial 
energy distribution of the mass-3 element. The stars in- 
dicate that the quantities are not independent from each 
other but appear under integral signs. 

Pai^a 2 with A\ ^ A 2 is found by integrating the non- 
thermal interaction rate over the time that takes A\ to 
slow down below the A 2 formation threshold E^ : 



Pa x ^,a 2 



dE Al 



^HeC-A!- 



..{E Al ) 



\dE Al /dx\ 



(6) 



Here, 7i4 He is the number density of ambient alpha parti- 
cles and o~a 1 ^a 2 is the cross section for the A 2 formation 
via A\ (a, N)A% with N — n or p. The energy loss rate 
dEA 1 1 dx of A\ is of great importance as it regulates the 
overall efficiency of the 9 Be producing chain; E A ^™ is the 
energy of A 2 once it is thermalized. For the tempera- 
tures that are of most interest to us, T < 20 keV, the 
thcrmalization rate is dominated by Coulomb stopping 
on electrons and positrons in the plasma. 

Equation ^ is to be averaged over the A\ energy dis- 
tribution /a 1 (E 1 2 1 ) = cr^ daAi/dE™ ; oa x is the cross 
section for producing A\ secondaries from He. In gen- 
eral, j Ax also depends on the energy of the primary which 
gives rise to A\. For the spallation products T and 3 He 
produced in N(a, X)T/ 3 He reactions, fs is largely inde- 
pendent of the incident nucleon energy, and we shall use 
the fit to fz,(E) as the sum of exponential factors 2l|; see 
also 



When it comes to A\ =6, 



3 Li and 6 He, we 



(5) 



assume a flat distribution of the differential cross section 
with respect to the CM scattering angle. 

When computing (Pg-yg)^ one should also account 
for the fact that 6 He may eventually decay (t6 Hc ~ 0.8 s) 
before interacting hadronically. This, however, has al- 
most no influence on ([5]). We find that even for temper- 
atures T as low as 1 keV more than 90% of 6 He survives 
the initial degradation from E™ to a kinetic energy of 
~ 1 MeV below which the cross section for ([2]) becomes 
too small. Turning to the last factor in (O, the proba- 
bility that a newly created 9 Be nucleus actually survives, 
we remark that in our BBN code the proton burning of 
thermalized 9 Be (as well as of 6 Li and 6 He) are already 
accounted for. Therefore, we identify (Pg^g) in ([3]) with 
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FIG. 2. Conversion probabilities Prie as a function of pho- 
ton temperature T for the various channels depicted; thermal 
destruction of 9 Be is not included. 



the "in-flight" destruction probability prior to thermal- 
ization. Finding (Pg^g) in an obvious modification of ([5]) 
we conclude that the survival probability is always well 
above 90% so that this factor is of no importance. On 
similar grounds, we can neglect any potential "in-flight" 
p-destruction of A = 6 secondaries. 

Figure [5] shows Pee for the various channels as labeled. 
The exponential drop in the number of electron-positron 
pairs leads to the sharp rise of Pee with temperature, 
with the peak at T ~ 20 keV, followed by a slow decline 
caused by the logarithmic growth of the stopping power 
at lower temperatures. 

We believe that our calculation of Pee is accurate to 
within a factor of two or three as it is largely based 
on experimental data. When available, we made use 
of the (inverse) cross section data collected in (23[. For 
the two nuclear reactions of most interest, 6 He(a 1 n) 9 Be 



25] 



and T(a,p) 6 He, we used the measurements of |2 
and [26| , respectively. The former reaction shows a large 
cross section above 100 mbn over a wide range of ener- 
gies, reaching a peak of 400 mbn at Eq ~ 5 MeV. The 
latter yields a maximum of 9 mb at E$ ~ 20 MeV and was 
approximated by a constant S- factor of S ~ 520 mb MeV 
for its inverse reaction 6 He(p, a)T. Previous investiga- 
tions of secondary 6 Li production have not accounted for 
the 6 He channel. This is well justified, because the rate 
for T(a,p) 6 He reaction is considerably less than that of 
T(a,n) 6 Li; the cross section for the latter as well as for 
3 He(a,p) 6 Li we take from [27j |. 

Constraints Having obtained Pb c , we now compute 
the actual output of 9 Be in the spirit of previous works 
on non-equilibrium BBN. The induced non-thermal con- 
tribution to the Boltzmann equation is of the form 
dn,9 Be /dt\ n . t h. = ^Hc^-s-s * Pbo, where r 4 ^ 3 is the de- 
struction rate of 4 He to A = 3 final states. Let X decay 
with some hadronic branching fraction B^ and denote by 
Ni(e) the number of energetic nucleons per V-decay with 




tx /sec 

FIG. 3. Contours of constant 9 Be (solid lines) and 6 Li (dashed 
lines) as a function of tx and e had Yx; B h — 1 and mx = 
1 TcV is assumed. 



kinetic energy e. Then, the spallation rate is given by 

r € had/2 



i 1 1 „ 



deN z (e)ai% 3 (e)v t (e) (7) 



where Chad is the injected hadronic energy per V-decay 
and -E^th ~ 20 MeV is the threshold energy in the A = 3 
producing reaction p/n + 4 He with cross section 174^.3. 
We do not need to account for other hadronic shower 
products such as long-lived mesons ir^, or Kl- For 
cosmic times t > 100s, they decay before interacting 
hadronically — unless the V- abundance is in excess of 



baryons [28[. We also find that photodissociation chan- 
nels such as 4 He(7, V)T/ 3 He and 7 Li( 7 , V) 6 He/ 6 Li are 
not important for the production of 9 Be. 

The time evolution of iVj is governed by the propaga- 
tion equation 



22] 



d t Ni = Ji- T hSC N t - d e (eT i>stop Ni) 



(8) 



where T^stop = —e~ 1 dEi/dt is the continuous stopping 
rate, Ti^ sc is the total hadronic scattering rate and Ji 
is the sum of all source terms. Up to secondary contri- 
butions from elastic and inelastic down-scattering pro- 
cesses, Ji is given by Jj ~ B^Q^e) /tx, where Qj(e) is 
the primary nucleon spectrum in the decay of V and Tx 
is the V-lifetime. Annihilating V are treated in complete 
analogy with the replacement r^ 1 — > nx((T aml f) where 
(cannv) is the annihilation cross section to hadrons. Fol- 
lowing 29|, 3(| we simulate a generic hadronic V-decay 
(and obtain Qi) by using the PYTHIA [31] model for 
e + e~ annihilations to hadrons. 

From inspection of TableUwe expect a C(l)-surviving 
fraction of non-thermally fused 6 Li or 9 Be only for 
T < 20 keV. For these temperatures, neutrons are not 
stopped before interacting hadronically. On the other 
hand, for protons with e < 1 GeV it is a reasonable ap- 
proximation to assume that the energy degradation term 



4 



in (J5J dominates. Then, employing a quasi-static equi- 
librium condition dtNi ~ to (JSJ (a good approximation 



for T in < r sc , r stop ) one finds 



eN p (e) 



r 



p,stop 



de' Q P (e'), 



eN n {e) eQ n (e) 



T ■ 



r 

1 n,s 



(9) 



Using © to obtain T 4 ^ 3 in ©, we solve the BBN net- 
work of Boltzmann equations for various combinations 
of tx and X-abundance Yx ■ The latter is normalized on 
the number density of baryons, Yx = nx/rib- 

Figure [3] shows contour lines of constant 9 Be/H and 
6 Li/H as a function of t x and Chad Yx- We have assumed 
B h = 1 so that e ha d = m x = 1 TeV. The (solid) 9 Be 
contours, when compared with the ones for 6 Li (dashed), 
illustrate that 9 Be is more sensitive to energy injection 
at earlier times, reaching its maximal abundance for 
Tx — 1.5 x 10 4 s. For larger lifetimes, as the primor- 
dial gas keeps rarefying, it becomes increasingly difficult 
to produce 9 Be from A = 6. For tx > 10 6 s the hadronic 
process looses further ground as neutrons decay before 
interacting strongly. At optimal tx, the current limits 
on primordial 9 Be probe the energy injection at a level 
of 10 MeV per nucleon, or 10 -12 GeV per en tropy . More 
elaborate treatments along the lines of 22, 29[ 3(| will 



certainly help to improve on the estimate of produced 
9 Be per unit injected energy. It is important to note, 
however, that the ratio 9 Be/ 6 Li, and thus their relative 
importance as constraints will not be affected by that. 
Another (competitive) constraint on hadronic energy in- 
jection is due to overproduction of D/H; for a comprehen- 
sive survey on general energy injection confer the works 
cited above. 

Can a chain of reactions similar to (JTJ lead to the fusion 
of 10 B? One way to increment A within the chain (JTJ by 
one unit is to replace one of its photonless recoil reactions 
with a radiative a-capture reaction, such as 6 Li(a,7) 10 B 
or 3 He(a,7) 7 Be followed by 7 Be(a,p) 10 B. The associ- 
ated cross sections are, however, only in the fib range, 
resulting in a 10 B conversion probability that is typically 
three orders of magnitude below the one of 9 Be. A more 
prospective avenue is to start off a sequence of transfor- 
mations with energetic 4 He rather than with T or 3 He. 
The former can be obtained in interesting quantities via 
elastic scatterings of shower-nucleons on ambient a par- 
ticles. This gives way to the formation of energetic A = 7 
from which 10 B (and 10 Be) can be fused. By explicitly 
modeling this path we find, however, that the A = 4 to 
A = 10 conversion probability Pb typically remains be- 
low 10 -10 so that 10 B is far from reaching observationally 
interesting levels. 

Conclusions In this letter we show that the produc- 
tion of "secondary" 6 Li [111 ] in non-equilibrium nucleo- 
synthesis scenarios can be followed up by a "tertiary" 
formation of 9 Bc/H in observationally relevant quanti- 
ties > 10~ 14 . The ensuing constraints on models of new 
physics are similar to those imposed by 6 Li, but do not 



suffer from the uncertainty associated with stellar deple- 
tion. Future observations of stars at lowest metallicities 
may strengthen this constraint even further, which sets 
9 Be apart from other light elements where such improve- 
ments appear unlikely 
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